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We study theoretically phonon-assisted spin relaxation of an electron confined in an elliptical quan- 
tum dot (QD) subjected to a tilted magnetic field. In the presence of both Rashba and Dresselhaus 
spin-orbit terms, the relaxation rate is anisotropic with respect to the in-plane field orientation. 
This anisotropy originates from the interference, at nonzero tilt angle, between the two spin-orbit 
terms. We show that, in a narrow range of magnetic field orientations, the relaxation rate exhibits 
anomalous sensitivity to variations of the QD parameters. In this range, the relative change in 
the relaxation rate with in-plane field orientation is determined solely by the spin-orbit coupling 
strengths and by the dot geometry. This allows simultaneous determination of both Rashba and 
Dresselhaus coupling parameters and the dot ellipticity from analysis of the angular dependence of 
the relaxation rate. 



I. INTRODUCTION 

Spin relaxation in semiconductor quantum dots (QDs) 
has recently attracted intense interest because of the 
possible use of the electron spin as a qubitji Quan- 
tization of orbital states in a QD due to the zero- 
dimensional (OD) confinement leads to the suppression of 
traditional spin-relaxation mechanisms (e.g., D'yakonov- 
Perel) that are dominant in continuous systems. Indeed, 
recent experiments on GaAs QDs in a magnetic field B 
have revealed extremely long spin-relaxation times (up 
to Ti = 170 ms at B ~ 1.75 T)^^^^d^>^^^ For mod- 
erate and high fields {B > IT), spin relaxation in 
QDs is dominated by phonon-assisted electronic transi- 
tions between Zeeman-split levels due to spin-orbit (SO) 
couplingi^ d^i^^d^d'^d^d^i^^ There are two distinct types 
of SO couplings, one originating from bulk inversion 
asymmetry (Dresselhaus coupling) and the other from 
structural inversion asymmetry along the growth direc- 
tion (Rashba coupling), that cause the admixture of or- 
bital states with opposite spins. "'^^ In the case of circular 
QDs in a perpendicular magnetic field, the Rashba and 
Dresselhaus terms mix different pairs of levels, and can, 
in principle, be distinguished if one such pair provides 
the dominant relaxation channel. This is the case, for 
example, when adjacent orbital levels, coupled via the 
Rashba term, are brought into resonance with changing 
magnetic field. However, in more realistic situations, 
deformations of the QD shape strongly alter the elec- 
tronic spectrunii^ and, in general, the effects of the two 
SO contributions are not separableJ^ 

An important distinction between the Rashba and 
Dresselhaus terms is their different symmetry proper- 
ties. The former, described by the Hamiltonian Hfj = 
otRi^xT^y — cTyTTa;), possesscs an in-plane rotational sym- 
metry, while the latter, Hd = aD{<^x''^x — Cj,7rj,), does 
not.^^- Here tt — — zV + eA, A being the vector poten- 
tial, <T is the Pauli matrix vector, and {old) is Rashba 
(Dresselhaus) coupling constant (we set Ti = 1). This 
lack of rotational invariance for Ho leads to an in-plane 
momentum azimuthal anisotropy in the presence of both 
SO terms that was recently reported in transport exper- 



iments in quantum wells. ^"-^^ In a magnetic field, the 
anisotropy arises due to the interference between Rashba 
and Dresselhaus terms in the matrix elements^^ in the 
presence of an in-plane field component. In QDs, this 
anisotropy reveals itself as a modulation of the spin re- 
laxation rate for different orientations of the in-plane 
field^^^i^^ 

Here we study the spin relaxation between Zeeman- 
split levels in elliptical QDs in a tilted magnetic field. We 
demonstrate that the interplay between SO interactions 
and QD geometry leads to dramatic changes in the relax- 
ation rate in a certain range of field orientations for which 
Rashba and Dresselhaus contributions undergo destruc- 
tive interference. Furthermore, in the vicinity of level an- 
ticrossings (see Fig. [T|), the SO contribution to the relax- 
ation rate factors out from the phonon one. This allows 
simultaneous determination of the parameters for both 
SO interactions and QD geometry from the azimuthal 
anisotropy of the differential (with respect to angle) re- 
laxation rate. 

The paper is organized as follows. In Section [TT] we de- 
rive electronic spectrum of elliptical QD in a tilted mag- 
netic field with both Rashba and Dresselhaus SO terms 
included. In Section [IIII the spin relaxation rate between 
lowest levels is evaluated. Numerical results for GaAs 
dots are presented in Section ITVl 



II. SPIN-ORBIT COUPLING AND ELECTRON 
STATES IN ELLIPTICAL QD IN TILTED 
MAGNETIC FIELD 

We start with the electronic spectrum in an elliptical 
QD in the presence of SO interactions subjected to a 
tilted magnetic field 

B = Bj^ + = B(^x sin 9 cos f + y sin 9 sin (p + z cos 9^ , 

(1) 

where 9 and ip are the tilt and azimuthal angles, re- 
spectively. The system is described by the Hamiltonian 
H = Hq + Hso + Hz, where Ho = ■n'^ /2m + Vc is the 
Hamiltonian of a 2D electron confined by the parabolic 
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FIG. 1: (Color online) (a) Energy levels and (b) relaxation 
rates ri2 and Fia for circular (dashed line) and elliptical 
(solid line) QDs vs. in-plane field component _B||. Arrows 
in panel (a) indicate relevant transitions, and circles indicate 
resonance regions. 



potential 



(2) 



u)x and LOy being the frequencies of the QD in-plane po- 
tential, Hso = Hr + Hd is the SO term, and Hz = 
^g* ^B<^ • B is the Zeeman term (m, g* and /is stand 
for the electron effective mass, g factor and Bohr mag- 
neton). For sufhciently strong 2D confinement, the or- 
bital effect of the in-plane field can be neglected and Hq 
depends only on the perpendicular field component via 
A — ■^(— J/, a;,0) (in symmetrical gauge). Using the 
transformation 



TTj; = aiPl - a2P2, TTy 

X = aiQi - 02(32, y = 



-mujibiQi + mLj2b2Q2 



bi 



■Pi 



-P2, 



(3) 



muJi muj2 

Pj and Qj {j = 1,2) being canonical momenta and 
coordinates, the Hamiltonian Hf) can be brought to 
the canonical form of two uncoupled oscillators with 
frequencies^ 



1 



(4) 

where ojc = eB±/m is the cyclotron frequency and the 
coefficients aj and bj are given by 



(5) 



with Dj = [(cj^ — -l-w^ijj^] ^ ^. In the absence of SO 
coupling, the spectrum represents two ladders of equidis- 
tant levels (for each spin projection) with energies 

<l.± = (^^1 + 1/2)^1 + ("2 + l/2)a>2 ± ^z, (6) 



where = \g* \i,BB is the Zeeman energy. 

The SO term Hso — Hd + Hr causes an admixture 
of the oscillator states \nin2)± with different orbital [rij) 
and spin (±) quantum numbers. In a titled field, the 
calculation of the SO matrix elements is carried out in 
two steps (see Appendix First, the spin operators 
a in Hso are rotated in spin space to align the spin- 
quantization axis along the total field B. Second, the or- 
bital operators tt in Hso are expressed via new canonical 
variables Pj , Qj . The expressions for the general matrix 
elements are provided in Appendix [XJ for the lowest ad- 
jacent levels with opposite spins corresponding to the L02 
ladder, they have the form 

±{Ol\Hso\00)^ 



hV2 



(62 COS 6 ± 02) cosif + i{a2 cos 9 ± 62) sin 93 
1^02 cos 6* ^ 62) cosiy9 — i(b2 cos6 ^ 02) sinip , (7) 



where Ij = {mujj)~^^'^ . These matrix elements explicitly 
depend on the tilt, 9, and azimuthal, ip, angles as well as 
on the QD geometry encoded in the coefficients aj , bj . In 
general, the magnitude of SO coupling is small compared 
to the level separation, a/lj <C UJj, and, accordingly, the 
SO-induced level admixture is weak. However, the level 
mixing gets strongly enhanced near the resonance, i.e., 
when the spacing between adjacent levels is of the order 
of the SO energy: lj2 — t^z ^ CiR,D/h (see Fig. [T]). This 
can be achieved, e.g., by varying the Zeeman energy with 
the tilt angle 9. The corresponding anticrossing gap. 



A-2| + (0l|ffso|00)_|, 
is evaluated from Eq. ([7]) as 



(8) 



'2 

a 



(62 cos^? -I- a2)^cos^ (y9 + (&2 + 0'2 cos6')^sin^ ip 



D 



+ 



(^2 cos 9 ~~ 02) sin^ + (62 ~ 0,2 cos 0) cos^ 

OiROiD I 2 , i2 



'2 



[cq + 62 ) sin 9 sin 2kp. 



(9) 



The gap magnitude is governed by the angle-dependent 
interference between Rashba and Dresselhaus terms. Im- 
portantly, in elliptical QDs, such interference depends on 
the dot geometry via the coefficients aj , bj . 

Thus, near the resonance, UJ2 — toz ^ A, the energies 
of the lowest excited states, 

£^2,3 = ^2 T ^ ^y{uj2-ujz)^ + A^, (10) 

acquire a strong angular dependence. At the same time, 
the phonon-assisted transitions between these states |j) 
{j — 2,3) and the ground state |1) are enhanced due 
to the strong admixture of constituent orbital levels 
(see Fig. [T]). As a result, the spin relaxation becomes 
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anisotropic with respect to the in-plane field orientation 
Lp. As we show below, the relaxation exhibits anomalous 
sensitivity of to the system parameters in a narrow range 
of ip where the SO terms interfere destructively. 



is the ratio of SO and orbital energies, while 77 = 
arg(+(01 |ifso |00)_) is the phase of the SO matrix el- 
ement. The parameter (3 characterizes the proximity to 
the resonance: 



III. PHONON-ASSISTED SPIN RELAXATION 

The transition rate between state |j) and the ground 
state |1) is given by 



QA 



£;,+^qa), (11) 



where the sum runs over acoustic phonon modes A with 
dispersion ojqx = c\Q, c\ being the sound velocity, and 
3D momenta Q — (q, (7z)- The transition matrix element 
is a product of phonon and electron contributions. 



(l|t/A|j)-M,(Q)(l|e^^^b-) 
where the phonon part, 

Ma(Q) = AA(Q)+^SA(Q), 



(12) 



(13) 



includes piezoelectric, A^, and deformation, Sa, 
contributionsj21 In the numerical calculations below, we 
include both longitudinal and two transverse piezoelec- 
tric acoustical modes. The details can be found in Ap- 
pendix |B] The electron matrix element can, in turn, be 
decomposed into a product of transverse and in-plane 
contributions, 



die' 



i) = fziqz)fijiq)- 



(14) 



The transverse contribution fz{qz) is determined by the 
2D confinement, assumed parabolic below, while the in- 
plane contribution /ij(q) is evaluated as follows. 

Not too far from the resonance region, I "''7'^^ I <^ 1 
it is sufficient to restrict calculations to the lowest four 
levels of the bj2 ladder. The states |2) and |3), corre- 
sponding to the anticrossing levels, are superpositions of 
unperturbed states, 



^00- 



|00)_ 



-^014 



101)- 



(15) 



while the ground state acquires a small admixture from 
the upper orbital of opposite spin, 



|l)=4oV|00> + 
The coefhcients d are obtained as 



4V-|oi)- 



(16) 



"01+ ~ "00- ^ 



1 



,(3) _ _ ,(2)* _ «^Ci 
"00- ^ "01+ ^ ^ \ ^ 



"00-1- — 



where 



-(Ol|i?so|00) 
0J2 + 



12(UJ2 + l-^z) 



(18) 



sinh/3 ~ (ti>2 — ^^z)/A 



(19) 



is the detuning in units of the anticrossing gap. In 
deriving Eqs. (|17p . we neglected higher-order terms in 
the SO coupling suppressed by the additional factor 
1^1 <C 1. In the resonance region /3 ~ 1, the coefficients 
c?(2) (^(3) ^ \^ j ^i^g levels 2 and 3 are well hybridized. 
Away from the resonance, A <C \u}2 — '^z\ ^ {1^2 + ^z), 
corresponding to \(3\ ^ 1, the eigenstates 2 and 3 are 
mainly determined by the unperturbed states |00)- and 
|01)+ (and vice versa) to the left (right) from the reso- 
nance region. 

In terms of the coefficients d, the 2D matrix clement 
takes the form 



/i.(q) = 4V-4^:/oooi(q) + <oV<-:/oioo(q), (20) 

with fm'n'mn = (to'ti' je^'^'" |mrt) . The relevant elements 
are (see Appendix [Cjl 



(1) M* 



0001 - Joioo - —f^^ 
with ( 



arg(q). We then obtain 



'(— a2 cos (l)+ib2 sin 0), 
(21) 



V2 



-10 



where 

^1^'^ = \ [(«2 + h2)d^iAol + («2 - 62)4^_4V: 



(22) 



1 



{a2 " 62)4iV4o+ + (02 + ^2)4-4^-1 • (23) 



Substituting Eq. 
ment. 



into the full transition matrix ele- 



(1|C/a|j) =MA(Q)/,(g,)/i,(q), 



(24) 



the scattering rate Tij is obtained by summing up over 
phonon modes in a standard manner-SiiSiiiiiiSiii^ (see Ap- 
pendix [BJ. The result is the product 



Ti, =Gi,(^)iy„ J = 2, 3 
where the geometric factor 



G 



0)1 2 



0)1 2 



(25) 



(26) 



is determined only by the SO-induced admixture of elec- 
tronic states, encoded in the coefficients d, while the 
phonon factor Wj describes the probability of phonon- 
assisted transitions between levels separated by energy 
Ej — El (see Appendix [B| . Note that Wj are nearly in- 
dependent of the SO coupling; in the resonance region. 



4 



the SO contribution to Wj is '-^ ^ 1 and can be ne- 
glected. 

Thus, near the resonance, the dependence of scattering 
rate T on the SO parameters and, accordingly, on the az- 
imuthal angle (/3, conies only through the geometric factor 
Gij. Remarkably, this factor can be extracted directly 
from the experimental data via the differential relaxation 
rate normalized to its value at some angle (e.g., ip = 0): 



AT, T,,{^) T,,{0) (^) - Gi, (0) 



r 



ri,(o) 



Gi,(0) 



(27) 



where the r.h.s. is independent of the phonon contribu- 
tion. 



IV. DISCUSSION AND NUMERICAL RESULTS 

Below we present our numerical results for spin re- 
laxation rates in a GaAs QD. Calculations were per- 
formed for a parabolic confining potential with cja, = 2 



meV and 



1.33 meV, corresponding to the elliptic- 



ity bj^l uiy — 3/2. For comparison, results for the circular 



QD (w^ 



2 meV) are presented too. We choose the 



parabolic transverse confinement with = 4 nm. The 
values of the SO constants, if not specified otherwise, in 
GaAs were taken as au. = 5 meVA and an = 16.25 
meVA, while the phonon parameters were taken from 
Ref. [li (see Appendix EI . 

In Fig. [T] we plot the lowest energy levels and relaxation 
rates for both circular and elliptical QDs as a function of 
the in-plane field B\\ at fixed = 6.0 T and <^ = 0. For 
chosen parameters, the level anticrossings at lo^ = ^z-, in- 
dicated by circled regions, for elliptical QDs are achieved 
at lower B\\ due to a weaker confinement along the y 
axis [see Fig.[lja)]. The relaxation rates and Fia are 
plotted in Fig. [T]Jb). The sharp increase in (Fis) is 
caused by a stronger SO-induced admixture of the |00)_ 

and |0l)_|^ states as one approaches to the resonance from 
the left (right). To the right of the resonance (wz > W2), 
is dominated by the orbital transition between states 
|00)^ and |0l)_|_; so is Fia to the left of the resonance 
(loz < W2). The flat B\\ dependence in these regions is 
because in a narrow 2D layer the orbital wave functions 
depend only on B^. Apart from the magnitude of F, the 
overall behavior is similar for circular and elliptical QDs. 

In Fig. [21 we plot the relaxation rate F12 as a function 
of in-plane field orientation Lp at several values of its mag- 
nitude in the resonance region. A strong azimuthal 
anisotropy is apparent: at certain angles, ri2 reaches 
minima that turn into maxima as B\ | sweeps through the 
resonance. For a circular QD, the extrema of ri2 occur 
at = ±7r /4 regardless of the values of a_R and a d [see 
Fig. mja)] . This anisotropy originates from the angular 
dependence of the anticrossing gap A((/9). Indeed, in the 
resonance region, \ijJ2 — ^z\/ A ~ 1, the SO contribution 




FIG. 2: (Color online) Angular dependence of the relaxation 
rate for (a) circular and (b) elliptical QDs for B±_ = 6.0 



T and several values of B\ 



to F12 takes the simple form 



Gl2((^) 




UJ2 - UOz 



2-1 1/2 



(28) 



For a circular QD, the expression 1 
11.22 



for the gap simplifies 



t^ 



A2 



2X]jX]i sin 2(p] 



(29) 



UJ2 

with \d = «£) sin^(6'/2)//2 and A_r = sin^(0/2)/Z2, 
so the extrema of A at (/? = ±7r/4 translate into the 
extrema of G. In contrast, in elliptical QDs, the angular 
dependence of F depends on the system parameters [see 
Fig. El^b)] . The additional asymmetry introduced by the 
QD ellipticity modifies the interference between Rashba 
and Dresselhaus terms and shifts the extrema away from 
±7r/4. The extrema of A, Eq. now depend on both 
the QD ellipticity and the SO parameters, and occur at 



tan 2(pe 



'■D "2 



hi 



(30) 



For the parameters of Fig. [H[b), — 106°. 

Such interplay between QD geometry and SO interac- 
tions suggests a way to simultaneously determine both SO 
couplings and QD ellipticity from the measured angular 
dependence of the differential relaxation rate AF/F, Eq. 
(I27p . In the resonance region, the phonon contribution 
Wj drops out, as does the prefactor in Eq. ([28)1 . so AF/F 
is determined solely by the anticrossing gap A(i^, ujx, ^y)- 
The angular dependence of AF/F is shown in Fig. [3] for 
several values of Rashba coupling a r, which can be var- 
ied, e.g., with an external electric field. For a circular 
QD, change in aR does not affect the minima positions 
aX if = — 7r/4, as discussed above; however, the modula- 
tion depth varies strongly [see Fig. [31^ a)]. Note that the 
ur dependence is nonmonotonic; the deepest minimum 
occurs for the almost complete destructive interference of 
the SO terms at aR/ao ~ \.&v?[0/2). 

Deviations from circular QD shape give rise to angu- 
lar dispersion of spin relaxation in the parameter space. 
For elliptical QDs, variations of ur shift the minima po- 
sitions of AF/F [see Fig. |31[b)]. In fact, the sensitiv- 
ity of AF/F to the system parameters is drastically en- 
hanced for if in the vicinity of the critical angles ^pe^ Eq. 
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This sensitivity arises from destructive interference be- 
tween the SO terms in a narrow domain of in-plane field 
orientations in the presence of asymmetric confinement. 
Note that such interplay between SO interactions and 
QD geometry cannot be captured by simplified descrip- 
tions of elliptical QDs using circular QDs with effective 
parametersi^ 



VI. ACKNOWLEDGEMENT 

This work was supported by the NSF under Grant 
No. DMR-0606509 and by the DoD under contract No. 
W912HZ-06-C-0057. 

APPENDIX A: SPIN-ORBIT MATRIX 
ELEMENTS IN TILTED FIELD 



FIG. 3; (Color online) Angular dependence of the differential 
relaxation rates for (a) circular and (b) elliptical QDs for sev- 
eral values of the coefficients an at B± = 6.0 T and = 13.2 
T (a) and B|| = 7.7 T (b). 

for which the destructive interference between the 
Rashba and Dresselhaus terms is strongest. Thus, a scan 
of the angular dependence of the experimental differen- 
tial relaxation rate in this narrow domain would enable 
an unambiguous extraction of both SO and QD geometry 
parameters. 

V. CONCLUSION 

In summary, we proposed a method for simultane- 
ous extraction of both the spin-orbit constants as well 
as the quantum dot shape from the angular anisotropy 
of the differential spin relaxation rate in a tilted field. 
The underlying mechanism is based upon the enhanced 
sensitivity of phonon-assisted spin-flip transitions to the 
system parameters in the vicinity of level anticrossings. 



Here we describe calculation of the SO matrix elements 
in elliptical QD in tilted magnetic field. It is convenient 
to work in the basis in which Zeeman term is diagonal. 
Therefore, we choose the spin quantization axis along the 
total field and perform the corresponding rotation of the 
Pauh matrices: 

fx (Jx cos9 cos(p — (Ty s'mLp + az sm9 cos(p, 
ay Ux cos 9 sin ip + ay cos tp + a^ sin 9 sin if, 
az —a^ sm9 + az cos9. (Al) 

The orbital variables are chosen to diagonalize the free 
Hamiltonian Hq with the help of transformation Eq. 
In this basis, the Hamiltonian Hq+Hz is diagonal in both 
orbital and spin spaces, with the eigenstates \ni,n2)± 
corresponding to two uncoupled oscillators with energies 
given by Eq. The calculation of matrix elements of 
Hso between these eigenstates is convenient to perform 
by utilizing rising = '^{j^ ~ ^'^^ lowering 

Cj = "^(^y^ operators. The corresponding non- 

zero matrix elements have the form 



±{ni,n2\Hso\ni,n2 - l)=p 



ri2f Or 

2 ^ h 



(62 COS0 ± a-i) cosiy9 -f z(a2 cos6' ± 62) sin 1^9 
02 cos 9 ± 62) cos -I- i (62 cos 9 =F a^) sin ip 



- — (62 cos (f + ia2 sm f) — t — — [a2 cos ip — 102 sm tp ) 
ni ( an 



smt 



± {ni,n2\Hso\ni — 1, W2)=f = a/ t ( T' "^^^ ^ T Qi) cos ip ± i(bi ± ai cos 6*) siniy9 

2 \ /i L 



+1 



(ai cos 6* q= bi) cos ip + i(— 61 cos6' ± ai) sinip 
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{ni,n2\Hso\ni - l,n2)± 



Or 
h 



(bi cos ip + iai sin ip^ 



.an , 

I— — ( fli cos If - 

h 



ibi sin ip] 



(A2) 



r 



Note that Hso does not couple states of different os- 
cillators. The second and fourth relations in Eqs. (|A2p 
describe SO-induced transitions without spin-flip that 
are absent in perpendicular field {6 = 0). Here, both 
Rashba and Dresselhaus terms contribute to the coupling 
between same levels. In contrast, for circular QD in per- 
pendicular field with LOx = — ^o, the two SO terms 
only couple different levels: 



with piezoelectric part containing LA and two transverse 
acoustical (TAI, TA2) modes. Since polarization direc- 
tions are given as 



e = X sm fq cos < 



-.TA2 



+ {ni,n2\Hso\ni,n2 - I) 
{ni,n2\Hso\ni,n2 - 1)+ = 



Or ln2UJ2 



n 



y sin 9q sin <j) + 
X cos 9q cos (f) + y cos 9q sin (; 
—X sin 4> + y cos (p 
[9q = arcsm(q/Q), (f> = arg(q)], one gets 



z cos 9q , 
) - z sin 6'q , 



(B5) 



(B6) 



.an n2Uj2 



where f2 = ^Jlo\^ uJ^jA: and aji^2 = ^1 ± ij-'c/2 (analogous 
relations hold for the ijJ\ ladder). 



where, as mentioned above, only LA mode is present for 
Sa. For Aa one obtains 

Ala(Q) = ^AoALA^sin2,/), 



AT^i(Q) = iAo^TA(2|-lj^/ 



2 sin 2(/), 



APPENDIX B: ELECTRON-PHONON MATRIX 
ELEMENTS 

1. Phonon matrix element Af(Q) 

The phonon-assisted transition rate is given by 



cos 20. 



(B7) 



r 



i,=27r^|(I|C/A|j)|2^(i?i-i?, 

QA 



In the above equations, Aq = e/114/K, and A\ — 

^12 

{2pVc\) . The sound velocities of the transverse 
and longitudinal modes are ctai = cta2 = cta and 
Cla, respectively. For GaAs, the parameter values are 
e/ii4 = 0.14 eV/A, Sq = 7 eV, cla ^ 5.14 x 10^ m/s, 
c\Q) (Bl) CTA = 3.03 X 10^ m/s, p = 5.31 g/cm^, and K = 12.79. 



where the operator of electron-phonon interaction U\ has 
the form^** 

C/a = {2pQcxVy'^^eAQ).+iQDQx)e''^'^{bQ^ + bl^^). 

(B2) 

Here, 



2. Expressions for Wj 



A 



QA 



{qxqy<^ 



QxqzGy 



qyQz 



(B3) 



is the amplitude of electric field created by phonon strain, 
and deformation potential Dqx contains only longitudi- 
nal acoustical (LA) component, Dq\ = Sx^la'^Oj with 
So being a constant of the deformation potential. Also, 
^QA (^Qa) creates (annihilates) phonon with dispersion 
'^QA = Qcx, V is the QD volume, R — (r, z) is 3D elec- 
tron radius vector, p is the crystal mass density, is the 
phonon polarization vector, cx is sound velocity, /114 is 
bulk phonon constant, k is the static dielectric constant. 
Accordingly, the phonon part Ma(Q) of the transition 
matrix element (l|J7A|j) includes both piezoelectric Aa 
and deformation Sa contributions: 



To obtain Eq. (P5|) . one has to perform integration over 
phonon modes in Eq. ljBl[) . The transition matrix ele- 
ment has the form {l\Ux\j) = Mx{Cl)fz{qz)fij{q), where 
A/a(Q) is described in the previous subsection, fz{qz) is 
provided in the next section, and fij is given by Eq. (|22p . 
Then, integration over eliminates (5-functions in Eq. 
(|B1[) . Subsequent </)-integration eliminates terms linear 
in F^ -''' . Remaining g-integration leads to Eq. (|25p where 
the phonon factor Wj is given by 



2tt 



xe 



I ^ f Ax{q)q^dq 




,(B8) 



with summation running over all phonon modes: defor- 
mation (Def), longitudinal and transverse acoustical: 



Ma(Q) = Aa(Q)+iSa(Q) 



(B4) 



Aoefiq) =^IAIa 



^LA 
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ALA{q) 
At A (q) 



6*^ 



1 CtA 



A CLA 



4 AK- 



K 


[Ai?,l 


' \ 1 
2 1 2 

- 9 j - q 


2-| 











3/2 



, (B9) 



where two transverse contributions are compacted as a 
single mode, AEj ^ Ej — Ei , and zero temperature has 
been assumed. 



APPENDIX C: ELECTRON MATRIX ELEMENTS 
IN ELLIPTICAL QD 

The form factor (g^) is calculated in the assump- 
tion that electron in the transverse direction is frozen 
on the lowest level. Its explicit expression for the trans- 
verse parabohc confinement Vz[z) = y^z^^ i^' fz{qz) = 



^1/2 



The explicit expression for the in-plane orbital matrix 
element, /„ 



im'n'mn = (m'n' e*'!'" wn) , is 



fm'n'mn — 



7'l„/l 



(0,1 cos (j) + ibi sin 0) 

v2 

— — ^ (02 cos (f> + ib2 sin c 



Li-H\^llW ),(C1) 



where L'^{x) is the associate Laguerre polynomial, and 
m! > m and n' > n. In deriving Eq. (|C1[) . the following 
identity has been used 

^,„'|g*'9[(aie-^'*+a2e'*)ct-H(Q2e-''*+aie'*)c]|^^ 

IT 



p— (ai+a2) I q /2 



xL™'-™[(ai-f a2)2;V]. (C2) 



^ D. Loss and D. P. DiVincenzo, Phys. Rev. A 57, 120 (1998) 
^ T. Fujisawa, D. G. Austing, Y. Tokura, Y. Hirayama, and 

S. Tarucha, Nature (London) 419, 278 (2002). 
^ R. Hanson, B. Witkamp, L. M. K. Vandersypen, L. H. 
Willems van Beveren, J. M. Elzerman, and L. P. Kouwen- 
hoven, Phys. Rev. Lett. 91, 196802 (2003). 

* J. M. Elzerman, R. Hanson, L. H. Willems van Beveren, B. 
Witkamp, L. M. K. Vandersypen and L. P. Kouwenhoven, 
Nature (London) 430, 431 (2004). 

^ M. Kroutvar, Y. Ducommun, D. Heiss, M. Bichler, D. 

Schuh, G. Abstreiter, and J. J. Finley, Nature (London) 

432, 81 (2004). 
® R. Hanson, L. H. Willems van Beveren, I. T. Vink, J. 

M. Elzerman, W. J. M. Naber, F. H. L. Koppens, L. P. 

Kouwenhoven, and L. M. K. Vandersypen, Phys. Rev. Lett. 

94, 196802 (2005). 
^ Y. Tokura, W. G. van der Wiel, T. Obata, and S. Tarucha, 

Phys. Rev. Lett. 96, 047202 (2006). 

* S. Amasha, K. MacLean, I. Radu, D. M. Zumbuhl, M. A. 
Kastner, M. P. Hanson, and A. C. Gossard, Arxiv eprint: 
[cond-mat/0607110 

^ X V. Khaetskii and Y. V. Nazarov, Phys. Rev. B 61, 12639 

(2000); 64, 125316 (2001). 
^° L. M. Woods, T. L. Reinecke, and Y. Lyanda-Geller, Phys. 
Rev. B 66, 161318(R) (2002). 

V. N. Golovach, A. Khaetskii, and D. Loss, Phys. Rev. 
Lett. 93, 016601 (2004). 
^2 D. V. Bulaev and D. Loss, Phys. Rev. B 71, 205324 (2005). 
C. F. Destefani and S. E. UUoa, Phys. Rev. B 72, 115326 



(2005). 

" V. I. Fal'ko, B. L. Altshuler, and O. Tsyplyatev, Phys. 
Rev. Lett. 95, 076603 (2005). 

J. Konemann, R. J. Haug, D. K. Maude, V. I. Fal'ko, and 
B. L. Altshuler, Phys. Rev. Lett. 94, 226404 (2005). 
P. Stano and J. Fabian, Phys. Rev. Lett. 96, 186602 (2006); 
Phys. Rev. B 74, 045320 (2006). 

See, e.g., R. Winkler, Spin-Orbit Coupling Effects in Two- 
Dimensional Electron and Hole Systems (Springer, Berlin, 
2003), and references therein. 

D. G. Austing, S. Sasaki, S. Tarucha, S. M. Reimann, M. 

Koskinen, M. Manninen, Phys. Rev. B 60, 11514 (1999). 

M. Valm-Rodrfguez, A. Puente, and L. Serra, Phys. Rev. 

B 69, 085306 (2004). 
2° J. B. Miller, D. M. Zumbuhl, C. M. Marcus, Y. B. Lyanda- 
Geller, D. Goldhaber-Gordon, K. Campman, and A. C. 

Gossard, Phys. Rev. Lett. 90, 076807 (2003). 
21 S. D. Ganichev, V. V. Bel'kov, L. E. Golub, E. L. Ivchenko, 

P. Schneider, S. Giglberger, J. Eroms, J. De Boeck, G. 

Borghs, W. Wegscheider, D. Weiss, and W. Prettl, Phys. 

Rev. Lett. 92, 256601 (2004). 

V. 1. Fal'ko, Phys. Rev. B 46, 4320 (1992). 

N. G. Galkin, V. A. Margulis, and A. V. Shorokhov, Phys. 

Rev. B 69, 113312 (2004). 

V. F. Gantmakher and Y. B. Levinson, Carrier Scatter- 
ing m Metals and Semiconductors (North-Holland, Ams- 
terdam, 1987). 



